The Origin of Nematic Order in the Iron-Based Superconductors 
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A theory of superconductivity in the iron-based materials requires an understanding of the phase 
diagram of the normal state. In these compounds, superconductivity emerges when stripe spin 
density wave (SDW) order is suppressed by doping, pressure or disorder, suggesting a possible role 
for magnetic fluctuations in the pairing mechanism. Experimentally, SDW order is often pre-empted 
by nematic order whose phase transition occurs at a slightly higher temperature but whose origin is 
yet to be resolved. One scenario is that the nematic order is unrelated to magnetism and reflects an 
orbital ordering of the iron 3<i-electrons. This ordering triggers the anisotropic magnetic coupling 
that produces stripe SDW order at a lower T. Another scenario is that stripe-SDW magnetism is the 
primary instability, but the magnetic transition occurs in two steps: first four- fold spin-rotational 
symmetry gets broken and the system develops spin- nematic order, and then time-reversal symmetry 
also gets broken, and full stripe-SDW order develops. Because spin-nematic order also induces 
orbital order, it is difficult to distinguish between the two scenarios experimentally. In this letter, 
we show that the magnetic scenario predicts that, at high doping, the system should undergo a 
transition at T < Tn into a new SDW phase in which four- fold rotational symmetry is restored. We 
have now observed a first-order transition to such a phase in Bao.76Nao.24Fe2As2, in good agreement 
with the spin-nematic model. 



INTRODUCTION 

There have been extensive investigations of the phase 
diagrams of the various iron arsenide and chalcogenide 
structures that display high temperature superconductiv- 
ity with critical temperatures up to 55 K [IHl]. In com- 
mon with other unconventional superconductors, such as 
the copper oxides, heavy fermions, and organic charge- 
transfer salts, superconductivity is induced by suppress- 
ing a magnetically ordered phase, which generates a high 
density of magnetic fluctuations that could theoretically 
bind the Cooper pairs. Whether this is responsible for 
the high transition temperatures has not been conclu- 
sively established, but it makes the origin of the magnetic 
interactions an important issue to be resolved [5] |6] . 

In nearly all the iron arsenides and chalcogenides, 
the iron atoms form a square planar net and the mag- 
netic order consists of ferromagnetic stripes along one 
iron-iron bond direction that are antiferromagnetically 
aligned along the orthogonal iron- iron bond [6j [7] . These 
systems are metallic and the Fermi surfaces, which are 
formed by the iron 3d electrons, are nearly cylindrical 
with hole pockets at the centre of the Brillouin zone and 
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electron pockets at the zone boundaries, all of similar 
size. In such an electronic structure, interactions be- 
tween electrons near the hole and electron pockets give 
rise to a spin density wave (SDW) order at the wavevec- 
tor connecting the two pockets [8] . This itinerant picture 
is consistent with the wavevector of the observed antifer- 
romagnetism, Angle Resolved Photoemission (ARPES) 
measurements of the electronic structure [9[ [10] , and the 
evolution of the dynamic magnetic susceptibility with 
carrier concentration [TTHT3] . 

However, any theory of the magnetic order also has to 
explain the structural transition which occurs at a tem- 
perature either just above or coincident with the SDW 
transition and lowers the symmetry from tetragonal (C4) 
to orthorhombic (C2). This is often referred to as nematic 
order, and the relation between nematicity, magnetic or- 
der, and superconductivity has become one of the central 
questions in the iron-based superconductors [TJJ [15] . 

At present, there are two scenarios for the development 
of the nematic order and its relation to SDW order. In 
the first, the structural order is unrelated to magnetism 
and is driven by orbital ordering as the primary insta- 
bility. The orbital ordering induces magnetic anisotropy 
and triggers the magnetic transition at a lower temper- 
ature by renormalizing the exchange constants [T6HT8] . 
This scenario is largely phenomenological, but there have 
been recent efforts to put it into a microscopic basis [19] . 

In the second scenario, the structural order is driven 
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by magnetic fluctuations, associated with the fact that 
striped SDW order can be along the x-axis (ordered mo- 
mentum is Qx = (0,7r)) or along the ?/-axis (ordered 
momentum is Qy = (7r, 0)). Theory predicts that the 
symmetry between the X and Y directions can be broken 
above the true SDW ordering temperature that breaks 
0(3) spin symmetry, i.e., the system distinguishes be- 
tween Qx and Qy without breaking time-reversal sym- 
metry [20]. The order parameter of this "Ising-spin- 
nematic state" couples linearly to the lattice, inducing 
both structural and orbital order. The magnetic scenario 
has been developed for itinerant [8[ [20] and localized [2TL 
[24] fermions, and the phase diagrams are rather similar 
in the two approaches. Below we use the fact that the 
systems we study are metals and use an itinerant ap- 
proach. 

It is difficult to distinguish the orbital and the mag- 
netic scenarios since many of the observable properties 
are identical. We have now identified both theoreti- 
cally and experimentally a new magnetic phase that can 
distinguish the two scenarios. By extending our pre- 
vious analysis of the itinerant model to higher carrier 
concentrations, we have discovered that the phase dia- 
gram is richer than previously anticipated, based on the 
Ginzburg-Landau analysis near T/v, and that there ex- 
ists, in some range of parameters, an additional phase at 
T < T/v, that restores C4 symmetry. The order parame- 
ter in this new phase combines Qx and Qy with equal 
weights. In the following, we summarize the results of 
these calculations, which are explained in more detail in 
the Supplementary Information, and show experimental 
evidence for the predicted phase in a hole-doped super- 
conductor, Bao.76Nao.24Fe2As2. Such a second transi- 
tion is highly unlikely in the orbital scenario because the 
breaking of the symmetry between Qx and Qy is the 
pre-condition for a magnetic transition. We view the ob- 
servation of the transition to an SDW state which does 
not break the symmetry between Qx and Qy as a strong 
indication that the nematic order is of magnetic origin 
rather than of orbital origin. 



THEORY 

The itinerant description of magnetism in iron-based 
superconductors is built upon the fact that the hole 
bands are centered around Qr = (0,0) and the electron 
bands are centered at Qx = (tt,0) and Qy = (0,7r), 
respectively (Fig. [lja)). A spin susceptibility is loga- 
rithmically enhanced at momenta connecting hole and 
electron pockets, and SDW order develops even if the in- 
teraction is weak. The SDW order parameter is in general 
a combination of the two vector components Ax and Ay 
with momenta (71", 0) and (0, 7r), respectively. For a model 
of perfect Fermi surface nesting (circular hole and elec- 
tron pockets of equal radii) and only electron-hole inter- 



actions, SDW order specifies |A^| 2 + I Ay 1 2 but not the 
relative magnitudes and directions of Ax and Ay . Away 
from perfect nesting, the ellipticity of the electron pock- 
ets and interactions between the electron bands break the 
degeneracy and lower the symmetry of the SDW order. 
Near T/v, an analysis within a Ginzburg-Landau expan- 
sion in powers of Aj and Ay shows that fourth-order 
terms select stripe magnetic order with either Ax ^ 0, 
Ay = 0, or Ay ^ 0, A x = [8] [20]. Such an order 
simultaneously reduces the lattice C4 symmetry down to 
C2. The order parameter in the stripe phase is shown 
schematically in FigjTJb). 

An issue that has not been discussed in detail until 
now is whether another magnetic ground state, in which 
both Ax and Ay are non-zero, may appear at a lower 
temperature, as a result of non-linear effects. This might 
happen either via a first-order transition, in which case 
the most likely outcome is the state in which | Ax | = | Ay | 
(see FigjTJc)), or via a second-order transition, in which 
case the second order parameter appears continuously 
and likely remains relatively small down to T = (see 
Fig0d)). 

To check for a potential second SDW transition, we 
extended the previous Ginzburg-Landau analysis to the 
entire temperature range by solving non-linear coupled 
mean-field equations for Ax and Ay and analyzing 
which solution minimizes the free energy. This has re- 
vealed new features in the phase diagram not previously 
identified. In particular, we find that SDW order with 
Ax = Ay, which breaks 0(3) spin symmetry but pre- 
serves lattice C4 symmetry, does emerge at low T, as the 
mismatch in hole and electron pocket sizes grows. 

We obtained this result by analyzing the minimal 
three-band model with one hole and two electron pock- 
ets. For simplicity, we considered parabolic dispersions 
with 
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where rrii are band masses, €0 is the offset energy, and \i 
is the chemical potential. 

The non-interacting Hamiltonian takes the form 



^0 = X^> k 4kcAka 



(4) 



where i — 1 — 3 label the bands, the summation over 
repeated spin indices a is assumed, and we shift the mo- 
menta of the fermions near the X and Y Fermi pockets 
by Q x and Qy, respectively, writing £x,k+Q x = 6,k, 
£r,k+Q Y = 6,k- 



3 




(b) $ $ $ (c) J^k J 1 (d) f $ f 
III ^SU^^i k 4 k 

A x =0, A Y ^0 |A X |=|A Y | * |A X |«^ Y |^0 

FIG. 1: (a) (upper panel) The band-structure with two cir- 
cular hole pockets at F and two electron pockets at X and 
Y, using the unfolded Brillouin zone with one Fe atom per 
unit cell. The arrows refer to two equivalent nesting wavevec- 
tors Qx = (tt, 0) and Qy = (0, tt) (lower panel) Possible 
magnetic ground states of the Fe-lattice: (b) the C2 antifer- 
romagnetic stripe phase with Ax = and Ay / 0; (c) a C4 
magnetic state, in which |Ax| = \ Ay \ 7^ 0, that is compatible 
with tetragonal lattice symmetry (other possible C± magnetic 
structures are discussed in the Supplementary Information); 
(d) magnetic order in which |Ax| << |Ay| / 0. Note that 
the Neel transition at low temperatures, T/v 2 , from phase (b) 
to phase (c) is first order, while the transition from phase 
(b) to phase (d) is of the second order in which an additional 
small component of Ax appears below Tn 2 • Our experiments 
are more compatible with scenario (c). 

The interaction term in the Hamiltonian Hint con- 
tains all symmetry- allowed interactions between low- 
energy fermions, which include inter- and intra-band 
scattering processes [25]. We present the explicit form 
of Hint m the Supplementary Information. The mean- 
field equations on Ax and Ay are obtained by in- 
troducing A x = (1/2JV) J2k 4,ka^£ c 2,k/3 and Ay = 
(1/2N) c\ kQ; <T a /3C3 5 k/3 and using them to decouple 
four-fermion terms into anomalous quadratic terms with 
inter-band "hopping", which depends on Ax and Ay. 
We diagonalized the quadratic form, re-expressed c^ka m 
terms of new operators and obtained a set of two coupled 
self-consistent equations for Ax and Ay. 

We solved the mean-field equations numerically as a 
function of two parameters, So and S2 (see Supplemen- 
tary Information for details). The parameter So = 2/i 
represents the mismatch in chemical potentials of the hole 
and electron pockets (So = when the electron and hole 
pockets are identical). S2 = Som(m x — m y ) / (2m x m y ) 
is proportional to the ellipticity of the electron pockets. 
We focused on the two SDW-ordered states- the antifer- 
romagnetic stripe state with Ax 7^ and Ay = 0, in 
which (^-symmetry is reduced to C2 and on the SDW 
state with Ax = Ay, in which (^-symmetry is pre- 
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FIG. 2: Phase diagram of Bai_ cc Na x Fe2As2. The blue points 
are the coincident antiferromagnetic and orthorhombic tran- 
sition temperatures, Tat, into the Ci phase, and the red point 
is the observed transition temperature into the tetragonal C4 
phase, all measured by neutron diffraction. The green points 
are the superconducting transition temperatures, T c , deter- 
mined from magnetization data. 

served. As we said, the two states are degenerate at zero 
ellipticity and perfect nesting, when S2 = So = 0. Once 
the ellipticity becomes non-zero, the stripe state wins 
immediately below the Neel temperature T/v. The C4- 
preserving state (with Ax = Ay) is a local maximum 
and is unstable at T < T/v. 

By solving the equations at lower temperature, we 
found that, at a finite So, the ^-preserving state also 
becomes locally stable below some T < T/v, and at an 
even lower T < T/v 2 , its free energy becomes smaller 
than that of the stripe phase, i.e., at T = T/y 2 the system 
undergoes a first-order phase transition in which lattice 
C4 symmetry gets restored, (see FigQc)). Because T/v 
falls as the Fermi surface mismatch £0 increases, the new 
C4-preserving phase in practice exists only in a narrow 
region of the phase diagram close to the suppression of 
SDW order. This behavior is illustrated in Fig. 2. We 
also analyzed a four-pocket model with two hole pockets 
and found another scenario for a second SDW transition. 
Namely, the AF stripe order Ay initially involves only 
fermions from a hole pocket which has higher density of 
states. Below some T < T/v, fermions near the remain- 
ing hole pocket and near the electron pocket at X, not 
involved in the initial stripe order, also produce SDW in- 
stability, and the system gradually develops the second 
order parameter |Ax|, which distorts the stripe AF or- 
der. The corresponding low T spin configuration is shown 
m Fig{TJd). In this case, however, the C4 symmetry re- 
mains broken at all temperatures. Our experimental data 
taken as a function of doping are more consistent with a 
first-order transition and restoration of C4 symmetry (see 
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below). It is possible, however, that the second scenario 
is realized under pressure [26] . 



EXPERIMENT 

We are currently exploring the phase diagram of 
Bai_ ;E Na :E Fe2As2 using neutron and x-ray powder 
diffraction, following our recent investigation of the 
potassium-doped compounds [1 j. Details of the synthe- 
sis and characterization of the polycrystalline samples 
and the powder diffraction measurements are given in the 
Methods section below. In both series, the addition of the 
alkali metal dopes holes into iron d-bands, which is equiv- 
alent to reducing the chemical potential, /i, in the model 
described above. This has the effect of reducing the tran- 
sition temperature into the stripe phase from 139 K, in 
the parent compound BaFe2As2, to at x ~ 0.25 in both 
the K-doped and Na-doped series. One unusual feature of 
both series is that the antiferromagnetic and orthorhom- 
bic transitions are coincident and first-order over the en- 
tire phase diagram [28], an observation that is quite un- 
ambiguous since both order parameters are determined 
from the same neutron powder diffraction measurement. 

We will publish a more detailed report on the whole 
phase diagram of Bai_ cc Na ;E Fe2As2 in a forthcoming pub- 
lication, but to understand its significance to this arti- 
cle, we show in Fig. 2 the transition temperatures into 
the antiferromagnetic/orthorhombic (AF/O) C2 (Tn) 
and superconducting phases (T c ) up to x = 0.3. In 
Bai_ ;E K ;E Fe2As2, it was possible to fit T/y vs x to a scal- 
ing form, i.e., to (1 — x/x c ) 2 ^ , with x c = 0.252 and 
/3 = 0.125 pQ. We do not have as many data points 
in the sodium-doped system, but the phase diagram is 
extremely similar so we have assumed that the same 
scaling curve applies here as well. This scaling has not 
been observed in the other iron-based superconductors 
and may reflect the reduced disorder in these hole-doped 
compounds, in which the doping occurs away from the 
iron arsenide layers. 

The only member of the sodium series to show a sig- 
nificant departure from the potassium compounds is at 
x = 0.24 close to the suppression of the AF/O order, 
which has a superconducting transition temperature of 
24 K. The transition into the conventional C2 phase oc- 
curs at 90 K. However, at T/v 2 = 39 K, there is a re- 
markable first-order phase transition, at which the sam- 
ple phase separates into regions with two different atomic 
and magnetic structures. Fig. 3(a) shows the orthorhom- 
bic splitting of the (112) nuclear Bragg peak at 90 K and 
its recombination at T/v 2 marking the onset of a reentrant 
tetragonal phase. Not all the sample transforms back to 
C4 symmetry. From a multi-phase Rietveld refinement, 
we estimate that approximately 40% of the sample re- 
mains in the orthorhombic phase below the transition 
(Fig. 4(a)). It is nevertheless clear that the whole sample 
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FIG. 3: First-order phase transition in Bao.76Nao.24Fe2As2. 
(a) A neutron diffractogram from the Wish Diffractometer 
showing the orthorhombic splitting of the (112) Bragg peak at 
Tn = 89 K and its recombination at the reentrant tetragonal 
phase at 38 K. (b) The lattice parameters, a (blue points) 
in the tetragonal phase of IA/mmm symmetry above Tn and 
below the reentrant transition, and a (red points) and b (black 
points) in the orthorhombic phase of Fmmm symmetry below 
Tn- The tetragonal lattice parameter is multiplied by y/2. 

is involved in the first-order transition even if only a frac- 
tion is transformed because, even within the orthorhom- 
bic phase, the in-plane lattice parameters relax below the 
transition (Fig. 3(b)). The in-plane lattice parameters 
show a reduction in the low-temperature C2 phase with 
an increase in the c-axis lattice parameter, i.e., the high- 
and low-temperature C2 phases are not identical. The 
C4 phase shows the opposite trend, which means that the 
in-plane lattice parameter is slightly higher than above 
Tn. The overall volume is roughly constant through the 
reentrant transition for both phases. 

Both the high-temperature and low-temperature C2 
phases have identical magnetic order with the moment 
in the plane, equivalent to the conventional AF stripe 
phase observed at other compositions. Within the high- 
temperature phase, the temperature dependence of the 
magnetic and orthorhombic order parameters are directly 



5 



(a) 1 



00 



80 

J 60 

o 

8 40 

CO 
Q_ 

20 




t 



************** *} 



- Tetragonal 
Orthorhombic 



J_ 



_L_ 



(b) 



10 20 30 40 50 
T(K) 



0.8 



if 
c 

CD 

E 
o 

2 0.4 

o 
"55 

I 0.2 
0.0 



I 


I 1 1 






I 


• — ■ — Orthorhombic 
— • — Tetragonal 1 
a-b/a+b 

I i I., 



0.002 



0.000 



30 60 
T(K) 



90 



FIG. 4: (a) The temperature dependence of the phase frac- 
tions of the tetragonal Ca phase (red points) and the or- 
thorhombic C2 phase (black points) below the reentrant tran- 
sition at 40 K derived from two phase Rietveld refinements, 
(b) The temperature dependence of the magnetic moment 
in the tetragonal C4 phase (red points), and the magnetic 
moment (black points) and orthorhombic order parameter 
5 — (a — b) I (a + b) (gold points) in the orthorhombic C2 
phase both above and below the reentrant phase transition. 



proportional (Fig. 4(b)), as we saw in the K-doped se- 
ries [28]. There is some evidence of phase competition 
above the reentrant transition since the moment in the 
low-temperature orthorhombic phase increases abruptly 
at the transition. The C4 phase is also magnetically 
ordered. An unambiguous refinement of the magnetic 
structure will require measurements on single crystals, 
but possible solutions that are consistent with the struc- 
tural space group (IA/mmm) and the lack of any in-plane 
magnetoelastic distortion are discussed in the Supple- 
mentary Information. They include the in-plane mag- 
netic order illustrated in Fig. 1(c), as well as variants, in 
which the moments are aligned along the c-axis. Since 
the theory does not include spin-orbit coupling, both so- 
lutions are theoretically allowed. 



DISCUSSION 

It is important to distinguish these new results from 
previous observations of a reentrant tetragonal phase in 
electron-doped compounds, such as BaFe2- x Co :c As2 [29] . 
All those transitions were within the superconducting 
phase and have been convincingly shown to result from 
the competition between superconductivity and antifer- 
romagnetic stripe order [30] [31]. The reentrant phase 
that we report here occurs well above the superconduct- 
ing transition temperature and so requires a different 
explanation. We are also unaware of any model of or- 
bital order that would predict a reentrant non-orbitally- 
ordered phase at lower temperature. Our results are 
broadly consistent with the itinerant model described 
above and are a natural consequence of magnetic order 
induced by the interaction between the hole and electron 
pockets suggested by the Fermi surface geometry. 

The relaxation of the lattice parameters below the 
reentrant transition implies that there is some phase 
competition between the C2 and C4 phases in the high- 
temperature phase that is relieved below T/v 2 • In princi- 
ple, this could be modeled theoretically by extending the 
mean-field theory calculations performed so far to the 
more microscopic theory. Such a theory would be neces- 
sary to study a C2 phase above the reentrant transition in 
which Ax 7^ Ay but both are non-zero. Nevertheless, 
the prediction that a C4 phase can become degenerate 
with the C2 phase only at higher doping when the hole 
and electron Fermi surfaces are not as well-matched in 
size, and that the stability of the C4 phase would be lim- 
ited to a very narrow region close to the suppression of 
antiferromagnetism is borne out by the new data. 

Our results provide strong evidence for the validity of 
an itinerant model of nematic order in the iron-based su- 
perconductors, in which the orbital reconstruction of the 
iron 3d states is a consequence of magnetic interactions 
induced by Fermi surface nesting. Whether nematic or- 
der, or at least strong nematic fluctuations, is a prerequi- 
site for superconductivity is another challenge to address 
in the future. 
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SUPPLEMENTARY INFORMATION 
Methods 

The synthesis of homogeneous single phase 
Bai_ ;E Na ;E Fe2As2 polycrystalline powders was performed 
using conditions that are similar to Bai_ ;E K :c Fe2As2. 
Further details of the synthesis can be found in Ref. pQ. 
Handling of all starting materials was carried out in a 
N2-filled glove box. Mixtures of Ba, Na, and FeAs were 
loaded in alumina tubes, sealed in niobium tubes under 
argon, and sealed again in quartz tubes under vacuum. 
The mixtures were heated at 850° C for one day and 
ground to homogeneous powders and then reheated for 
3 days at 850° C. This last step was performed twice. 
Initial characterization of the the final black powders 
was performed by x-ray diffraction, magnetization 
measurements, and Inductively Coupled Plasma (ICP) 
elemental analysis. Powder diffraction experiments were 
performed on three time-of-flight neutron diffractome- 
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ters, HRPD (x = 0, 0.24) and WISH (0.15, 0.24) at 
the ISIS Pulsed Neutron Source (Rutherford Appleton 
Laboratory) and HIPD (0.1, 0.2, 0.3) at the Lujan 
Neutron Scattering Center (Los Alamos National Lab- 
oratory). Synchrotron x-ray diffraction measurements 
were performed as a function of temperature for x = 
0, 0.1, 0.15, 0.2 at 11BM-B at the Advanced Photon 
Source of Argonne National Laboratory. The nuclear 
and magnetic structures and internal parameters were 
determined by the Rietveld technique using GSAS [2] 
and EXPGUI El. 



Theory of the Reentrant C4 Phase 

The magnetic phase in most parent compounds of 
the iron-based superconductors is the stripe spin-density 
wave order with momentum either Qx = (0, n) or Qy = 
(7r, 0) in the unfolded Brillouin zone, in which there is 
one iron atom per unit cell [4j |5] (see Fig. 1(b) in the 
main manuscript). These two wavevectors connect the 
hole pockets at the center of the Brillouin zone and the 
electron pockets at the zone boundary, along the two or- 
thogonal iron-iron bond directions. The stripe magnetic 
ordering breaks both 0(3) spin-rotational symmetry and 
C4 lattice rotational symmetry (by selecting either Qx 
or Qy), and is often preceded by a "nematic" phase in 
which C4 symmetry is broken, but 0(3) rotational sym- 
metry remains unbroken. 

In general, the geometry of iron pnictides allows more 
complex orders in which both Aj and Ay are present. 
Previous analysis by two of us [6 have shown that near 
T/v, when Ginzburg-Landau theory is valid, the system 
definitely prefers a stripe order in which only the order 
parameter with Qx or Qy is non-zero. Here we extend 
the previous analysis to lower T and study the magnetic 
order between T = Tn and T = 0. Our results show 
that the phase diagram of the system is more complex 
than previously thought. In particular, we find that at 
some range of dopings, the system undergoes a first-order 
transition, upon lowering T, into a phase in which there 
is a two-component order parameter with equal magni- 
tudes of the components with Qx and Qy, and the four- 
fold lattice rotational symmetry is restored (we label this 
phase as a O4 phase). The first order phase transition be- 
tween the O2 stripe and O4 phases is consistent with the 
observed by neutron diffraction experiments reported in 
the main text. 

We consider the minimal three-band model with the 
hole pocket T at the center of the Brillouin zone and two 
electron pockets X and Y at Qx and Qy, respectively. 
For simplicity, we consider parabolic dispersions with 



2m x 



k- 

2m v 



-/i, and 



£r,k = £0 - 2^ - (x.k+Qx = —Go - 

k 2 k 2 

£y,k+Q Y = -£o + 2^T y + V' where m, denotes the 

band masses, 6q is the offset energy, and \i is the chemical 



potential. Near the Fermi energy and for small elliptic- 
ity, the dispersions can be approximated by £r,k = — 
£x,k+Q x = f _ ^0 + h cos 2/9, £y,k+Q Y = f - - £2 cos 2(9, 
with £q = 2/i, 82 — £(yrn(m x — m y ) / {2m x m y ) 1 and 
= tan -1 k y /k x [9] . 

Electrons with spin a of the band i are created by the 
operators cj~ kQ; , and free-fermion part of the Hamiltonian 
has the form 



n 



i,k 



(5) 



Here the summation over repeated spin indices is as- 
sumed, and we shift the momenta of the fermions near 
the X and Y Fermi pockets by Qx and Qy, respectively, 
i.e., write £x,k+Q Y = fx,k, £r,k+Q Y = fy,k. 

To shorten presentation, we restrict the interacting 
part of Hamiltonian to the interaction in the spin channel 
with momenta near Qx and Qy, i.e., to 



j [/spin ^ ^ ^i,q • 



(6) 



where s^ q = c F k+qQ ,cr a/ #Ci 5 k/3 is the electronic spin 
operator, and a a p are Pauli matrices. The coupling [/ sp in 
is the sum of density-density and pair-hopping interac- 
tions between hole and electron states ([/spin = U\ + Us 
in the notation of Ref. [H]), where 

u i4",a c r,<xC X (3 c x ,(3 = -^c^a^cx,^ • Cx, 7 0V c r,<5 
+(■■■) 

TT t t ^3 + t 

U 3 c r,a C X,aCr,{3 C X,P = ~ c r> a c^ C X,/3 ■ C r?7 <T 7( 5Cx,(5 

+(•••) (7) 

and the dots stand for the terms with 5 aj p5 7j §, which 
only contribute to the CDW channel. The couplings U\ 
and Us do depend on the angle along the electron pock- 
ets [7], but for our purposes this dependence may be ne- 
glected, i.e., [/ S pin can be approximated by a constant. 
Once [/ sp in exceeds some critical value (which gets larger 
when So and 62 increase), the static magnetic susceptibil- 
ity diverges at (0,7r) and (7r,0), and the system develops 
long-range magnetic order. 

To understand what kind of magnetic order wins be- 
low T/v we introduce the two spin fields A( X ,y) = 

c r ka (7 ^C(x ) Y),k/3- We apply Hubbard- 
Stratonovich transformation, integrate out fermions, ob- 
tain the action S[Ax, Ay] in terms of Ax and Ay, 
use saddle-point approximation dS/dAi = 0, and solve 
a set of coupled saddle-point equations for Ax and Ay 
A straightforward way to perform this calculation is is to 
introduce the 6-dimensional Nambu operator 



k — ^ c r,kt 



r,k| 



X,kt 



4 



) (8) 
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Applying the Hubbard-Stratonovich transformation and 
evaluating the products of the Pauli matrices, we obtain 
the partition function in the form [SJ [TU] : 



(9) 



with the action 



5[*,Ai] = -/^-i^ fe+ 2 f ( 

Jk ^spin Jx 



A^ + Af.) (10) 



Here A$ = |Ai|, and the Green's function Q k 1 is given 
by: 



S k 1 = Go,l ~ v 

where the free-fermion term is 

Gr,k 
G x ,k 
G Y , k 

and the interaction term is 



(11) 



(12) 



(13) 



Here we introduced the 2x2 matrices G^k — G\^L and 
A^ = Ai -cr, where I is the identity matrix. The functions 
Gil ~ ^ ~ <^> k are tne non-interacting single-particle 
Green's functions for T, Y, and X fermions. 

It is now straightforward to integrate out the fermions, 
since the action is quadratic in them, and obtain the 
effective magnetic action: 





( ° 


-A x 


-Ay 


H 


i ~ Ax 










K~A Y 









5 eff [A x , A Y ] = -Tr In (1 - Q 0ik V)- 



x spm Jx 



(A^ + A y ) 
(14) 

Here Tr (• • •) refers to the sum over momentum, fre- 
quency and Nambu indices. In contrast to the previous 
studies [6j [10] , we do not perform a series expansion in 
powers of A 2 but analyze the full non-linear saddle-point 
(mean-field) solutions j^- — for the two cases: (i) a C2 
stripe phase, in which we set Ax ^ 0, Ay = and (ii) 
a C4 phase in which we set Ax = Ay = A For the case 
(i), the mean- field equation has the form 



1 = 2U S p in Ax ^2 



1 



A 2 _ ri-l (-i-l 



while for the case (ii) we have 
1 — 2£/ S p2 n A ^ ^ 



^-G-lG-) k + ^G-\G 



(15) 



(16) 



In both cases the sum over Matsubara frequencies can be 
evaluated exactly. For stripe magnetic order, the mean- 
field equation becomes 

1 = -2U spm A x £ (17) 





FIG. SI: (color online) Calculated magnetic phase diagram 
as a function of the ellipticity 62 at a finite mismatch between 
the electron and hole pockets, So = 0.01 (a) and as a function 
of the mismatch So for the finite ellipticity S2 — 0.01 (b). The 
parameter So increases with increasing doping. The squares 
denote the Neel temperature for the stripe phase. The trian- 
gles denote the temperature below which the C4 phase wins 
over the stripe phase, although the two phases remain nearly 
degenerate in energy. The solid and dashed curves are guides 
to the eye. We used U s dw — 0.8eV, eo = 0.2eV, and m — leV. 
In Fig. |S1| 3, the energies of the C± and the stripe phases are 
very close for all T, so Ca phase may actually win at low T 
in a more generic model. 



where 



1 



£i,2k = ^ ( Sr,k + £x,k+ Qx ± V (&,k - ^,k +Qx ) z + 4A X 



^ 2 +4A 2 
(18) 



. For the C4 phase, we obtain 

C^Hk 



1 



■ £x,k+Q Y ) / (^llk) 



(^llk — ^22k) (^llk — ^33k) 



(£?22k - £x,k+Q Y ) / (^22k) 
(^llk — ^22k) (^22k — ^33k) 
(^33k - £x,k+Q Y ) / (^33k) 



(^llk — ^33k) (^22k — ^33k)_ 



(19) 



where the energies En (i = 1—3) are the three solutions of 
the cubic equation A 2 (uj - £x,k+Q x ) + A 2 (cj - £y,k+Q Y ) - 
- Cr,k)(w - ^x,k+Q x )(^ - (y,k+Q Y ) = 
We solved these equations numerically together with 
the equation for the chemical potential, for different val- 
ues of the chemical potential mismatch So and ellipticity 
parameter 5 2 - Like in previous analysis [8], we find 
that the actions for C2 and C4 phases are degenerate 
at zero ellipticity and for equal sizes of the electron and 
hole pockets (So = S2 = 0) . Once ellipticity becomes non- 
zero, the C2 wins near the Neel temperature. Within the 
Ginzburg-Landau expansion to order A^ y , the lower 
free energy of the C2 phase is the consequence of the 
fact that the ellipticity generates the term C|Ax| 2 |Ay| 2 
with positive coefficient C, which increases the energy 
of the C4 phase but does not affect C2 phase. Going 
beyond Ginzburg-Landau approximation, we found that 



for small enough ellipticity, the solution of Eq. (19) re- 
emerges below some T < T/v and, below this T, C4 phase 
becomes a local minimum. Furthermore, in some range 
of A2, at even lower T < T/V2, the free energy of the 
C4 phase becomes slightly smaller than that of the C2 
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Model A 

X X M = (0,0,-M) 



o 


o 


o 


M = (0,0,M) 
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Fcmm'm' 
a' = ^2a 



Model B 



Pp4/mmm 
a' =\/2a 



At z + 1/2 each sign is reversed 



FIG. S2: The two magnetic models for the reentrant tetrago- 
nal phase that best fit the data. In Model A, the circles and 
crosses represent moments that are parallel and anti-parallel, 
respectively, to the c-axis. 



phase, i.e., the system undergoes a first-order transition 
from C2 to C4 magnetic phase, and lattice C4 symmetry 
gets restored. 

At Sq = the region of 82 where C4 phase wins is 
exponentially small. However, at a finite 5q? this region 
widens up. We show the results for a particular So in 
Fig. [Slja) and for a given ellipticity as a function of So 
in Fig. Sl/b). A non-zero So means that the sizes of 
electron and hole pockets are non-equal, i.e., that there 
is a finite amount of doping. The implication of this 



result is that, at a finite doping, as the temperature is 
lowered, the system first orders into a stripe C2 phase in 
which four-fold lattice rotational symmetry is broken (X 
and Y directions become non-equivalent), and then, at 
a lower T, it undergoes a first order transition into the 
C4 phase in which four-fold lattice rotational symmetry 
is restored. This is consistent with our experiment. 



Possible Magnetic Ordering in the C4 Phase 

High resolution neutron diffraction on HRPD (ISIS 
Facility) shows that the atomic structure of the C4 
phase has the same IA/mmm space group as the high- 
temperature paramagnetic phase, but coexists with the 
C2 phase with the orthorhombic Fmmm space group be- 
low T/v 2 = 39 K (see Fig. 3 of the manuscript). Since 
the ordered moments are so small (~ 0.4 to 0.8/ie), the 
magnetic peaks are too weak to be observed on HRPD, 
but are measurable on the Wish diffractometer (ISIS Fa- 
cility), where they are comparable in intensity to Bragg 
peaks due to Fe2As and other impurity phases. This 
makes a completely unambiguous refinement of the mag- 
netic structure, and even the direction of the magnetic 
moments, impossible. However, the symmetry imposed 
by the atomic structure and the constraints imposed by 
the data allow us to focus on the models discussed below. 

Most simple magnetic models can be ruled out due to 
their inability to fit the large increase in intensity in the 



(101)f peak that coincides with the C4 phase (the sub- 
script denotes that it is indexed based on the F-centered 
a' = V2a cell). Only two models fit the data satisfacto- 
rily (Fig. |S2[ ). In the first (A), the magnetic moment is 
oriented in the c direction, with stripes in either the a or 
b direction (which are equivalent). In this case the mag- 
netic space group is Fcmm'm' . In the second model (B), 
there are no stripes; instead the moments lie in the ab 
plane, spacegroup PpA/mmm (Pc^/mmm with respect 
to the IA/mmm cell). Moments adopt orientations in the 
(110) p directions, with the C4 axis mandating that each 
moment is orthogonal to its nearest neighbor (this is the 
model illustrated in Fig. 1(c) in the main manuscript). 
In both A and B models, each layer is coupled antiferro- 
magnetically to its neighbor, which is determined by the 
intensity of the magnetic (10/)f reflections. (In the AFM 
case (10/) = for I = 2n and in the FM case (10/) = 
for I = 2n + 1. We observe the former.) 

The fits to the two co-existing phases observed below 
39 K, i.e., the C4 phase (models A or B) and the C2 
stripe phase, are compared in Fig. [S3] Model A pro- 
duces a slightly better overall fit than model B. However, 
the weak signal-to-noise ratio and the overlap with peaks 
from the coexisting C2 magnetic phase as well as weak 
impurity phases preclude an unambiguous solution. The 
magnetic moments reported in the manuscript (see Fig. 
4 in the main manuscript) for the C4 and C2 phases are 
based on a free refinement of the moments for the C4 
model A and the standard C2 model. 



Relationship of Magnetic Ordering to the Theoretical 
Electronic Model 



The atomic structure of the C4 phase is tetragonal, 
with no evidence of orthorhombic strain even in the high 
resolution measurements on HRPD. This implies by sym- 
metry that J a = J5, which is consistent with Model B, 
whose magnetic space group is also tetragonal. However, 
the magnetic structure of Model A is still apparently or- 
thorhombic even though the moments are along the c- 
axis. This cannot be reconciled with an orbital order 




d(A) or (A) 

FIG. S3: Fits to the principal magnetic peaks at various tem- 
peratures: a) 106 K (tetragonal non-magnetic), b) 46 K (or- 
thorhombic phase only), c) 2 K (tetragonal model A + or- 
thorhombic), and d) 2 K (tetragonal model B + orthorhom- 
bic). All plots are on same scale, error bars are d-dependent 
in neutron time-of- flight measurements. 



10 





cos(nh) 






cos(nk) 




COS(TTh) + C0S(7Tfc) 


X 


XXX 


X \~ b 
1 


X 


O X o 


X 


XXX 











X 


0X0 


X 


00 


X 


XXX 


x + 


X 


O X o 


X 


= XXX 











X 


0X0 


X 


■O 


X 


XXX 


X 


X 


0X0 


X 


XXX 

|/W,| = \M\ + |M| = 2[M| 
\M 2 \ = \M\ - \M\ = 




cos(tt/i) 






cos(nk + q> 




cos(Trh) + cos(n7( + <p) 


X 


XXX 


x r 


x 


O X o 




-m,x x x x x 





O 





X 


O X o 


X 


M 2 o 


X 


XXX 


x + 




O X o 




= X x x x X 





O 


M 




0X0 




o o ° 


X 


XXX 


X 




O X 




X x X x x 






-M 




5M 


-5/W 


\M,\ = \M\ + \5M\ 
\M 2 \ = \M\ - \5M\ 



that the striping alternates directions between layers, so 



FIG. S4: Upper panel: A combination of the two equivalent 
arms in model A gives a checkerboard with half moments 
equal to zero. Lower panel: Adding the arms out of phase 
yields incomplete cancelation of moments in the SDW model. 



model, in which the stripes are a result of J a ^ J5, but it 
is also not consistent with the itinerant model, in which 
the a and b axes are equivalent. However, tetragonal 
symmetry can be restored if we assume that, instead of 
separate domains of a and 6-axis stripes, the magnetic 
structure involves a simultaneous modulation in both di- 
rections, i.e., a double-Zc structure. 

Writing the two forms as a spin-density waves with 
modulations cosfjrh) and cos(7rk) shows that adding the 
waves interfere such that half of the moments double and 



the others sum to zero (Fig. S4 upper panel). The nodes 
on every other Fe atom in a checkerboard fashion remove 
striping and recover tetragonal symmetry. In an itiner- 
ant magnetic model, a node implies that there is actually 
zero moment centered on the nodal atom (unlike localized 
moment systems where nodes in a magnetic model imply 
a fluctuating spin and therefore contribute a significant 
entropy to the system) . This is physically allowed but an 
alternative choice without nodes is to change the phasing 
between waves, e.g., cos(7rk) + cos(7rk + </>), which scales 
the amplitude of one direction because cos(irk + </>) = ±5 
where < S < 1. This gives a situation where \M\\ = 
\M\ + S and \M 2 \ = \M\ — S . This gives a checker- 



board of alternating large and small moments (Fig. S4 
lower panel). This can be modeled in tetragonal space 
groups without a proper C4 axis, such as PA2 and P4 
which allow for two inequivalent coplanar Fe atoms. We 
modeled this arrangement (model A2) in the spacegroup 
PpA'2\d (#114). This maintains the striping, but there 
is an intractable difference between models A and A2. 
The 4 axis preserve tetragonal symmetry by requiring 



that there is no preferred direction (left panel, Fig. S5). 
Thus, the only situation that is fully AFM coupled be- 
tween layers is when M2 = 0. This makes unobserved 
reflections (102n) non-zero. However, situations where 
Mi ^> M| (right panel, Fig. S5) predict intensities for 



these peaks that are on the same scale as the noise. Oth- 
erwise the two models are identical. In short, Model A 
is a more elegant solution to the magnetic data than A2. 
However, a final determination of the magnetic structure 
will have to wait for single crystal measurements. 
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FIG. S5: Left panel: Comparison of the differences between 
layers in the magnetic models A and A2 described in the text. 
Right panel: An example of a refinement of Model A2. 
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